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Abstract 

The aim of this paper is to show an existence theorem for a kinetic model of 
coagulation-fragmentation with initial data satisfying the natural physical bounds, 
and assumptions of finite number of particles and finite L p -norm. We use the no- 
tion of renormalized solutions introduced dy DiPerna and Lions in [3], because of the 
lack of a priori estimates. The proof is based on weak-compactness methods in L , 
allowed by L p -norms propagation. 



1 Introduction 

Coalescence and fragmentation are general phenomena which appear in dynamics 
of particles, in various fields (polymers chemistry, raindrops formation, aerosols, ...). 
We can describe them at different scales, which lead to different mathematical points 
of view. First, we can study the dynamics at the microscopic level, with a system of N 
particles which undergo successives mergers/break ups in a random way. We refer to 
the survey [1] for this stochastic approach. Another way to describe coalescence and 
fragmentation is to consider the statistical properties of the system, introducing the 
statistical distribution of particles f(t,m) of mass m > at time t > and studying 
its evolution in time. This approach is rather macroscopic. But we can put in an 
intermediate level, by considering a density / which depends on more variables, like 
position x or velocity v of particles, and this description is more precise. Here, we start 
by discussing models with density, from the original (with / = f(t,m)) to the kinetic 
one (with / = f(t,x,m,v)), which is the setting of this work. 

Depending on the physical context, the mass variable is discrete (polymers forma- 
tion) or continuous (raindrops formation). It leads to two sorts of mathematical models, 
with m 6 N* or m e]0, +oo[, but we focus on the continuous case. To understand the 
relationship between discrete and continuous equations, see [16j. 
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1.1 The original model 



The discrete equations of coagulation have been originally derived by Smoluchowski 
in [21], [22], by studying the Brownian motion of colloidal particles. It had been extended 
to the continuous setting by Miiller [20], giving the following mathematical model, 
called the Smoluchowski 's equation of coagulation: 

^(t,m)=Qt(f,f)-Q:(f,f), (t,m) G]0,+oo[ 2 . (1.1) 

This equation describes the evolution of the statistical mass distribution in time. At 
each time t > 0, the term Q^if, f) represents the gain of particles of mass m created 
by coalescence between smaller ones, by the reaction 



{m*} + {m — m*} — > {m}. 



The term Q c (/, /) is the depletion of particles of mass m because of coagulation with 
other ones, following the reaction 



Namely, we have 

Qt(fJ)(t,m) 

Q-{fJ)(t,m) 



{m} + {m*} — > {m + m*}. 



A(m*, m - m*)f(t, m*)f(t, m - m*)dm* 



A(m, m*)/(t, m)f(t, m*)dm* 



where Aim, m*) is the coefficient of coagulation between two particles, which governs 
the frequency of coagulations, according to the mass of clusters. In his original model, 
Smoluchowski derived the following expression for A: 



A(m, m* 



m V3 + m *l/3\ ( m -V3 + m *-l/3 



;i. 2 ) 



In many cases, coalescence is not the only mechanism governing the dynamics of 
particles, and other effects should be taken into account. A classical phenomenon 
which also occurs is the fragmentation of particles in two (or more) clusters, resulting 
from an internal dynamic (we do not deal here with fragmentation processes induced by 
particles collisions) . Th is binary fragmentation is modeled by linear additional reaction 
terms in equation fll.ip . namely 



Q+(f)(t,m) 



+oo 



B(m , m)f(t, m)dm! , 



i r m ' 

Q](f)(t,m) = -f(t,m)B 1 (m), where B x {w!) = / B(m', 

K ^ JO 



m)dm. 



The function B(m' } m) is the fragmentation kernel, it measures the frequency of the 
break-up of a mass m! in two clusters m and m — m, for m < m . So, at each time t, 
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the term Q~f(f) is the gain of particles of mass m, resulting from the following reaction 
of fragmentation: 

{m!} — > {in} + {m! — in}, 

whereas Qj(f) stands for the loss of particles of mass m, because of a break-up into 
two smaller pieces, by the following way: 

{m} — > {m*} + {m — m*}, with m* < m. 

Thus, the continuous coagulation-fragmentation equation writes 



-±(t, m) = Qt(f, f) - Q7(f, f) + Qj{f) - Q](f), (t, m) G]0, +oo[ 2 . (1.3) 

In the 90's, many existence and uniqueness results have been proved about this problem, 
see for instance [23], [24], or [13] for an approach by the semigroups of operators theory. 
These results are true under various growth hypotheses on kernels A and B, but these 
assumptions often allow unbounded kernels, which is important from a physical point 
of view. 

However, this coagulation-fragmentation model do not take the spatial distribution 
of particles into account. This leads to "spatially inhomogeneous" mathematical mod- 
els, where the density of particles f(t,x,m) depends also of a space variable x G M 3 . 



1.2 Spatially inhomogeneous models 

A first example consists of diffusive models, corresponding to the situation where 
particles follow a Brownian motion at the microscopic scale, with a positive and mass- 
dependent coefficient of diffusion d{m). From a physical point of view, it implies that 
particles are sufficiently small to undergo the interaction with the medium, i.e the 
shocks with the molecules of the fluid in which the particles evolve. In the statistical 
description, a spatial-laplacian appears, giving the diffusive coagulation-fragmentation 
equation: 

df 

Mt, x, m) - d(m)A x f(t, x, m) = Q+{f, f) - Q~ (f, f) + Q+(f) - Q] (f), 

dt (1.4) 

(t,x,m) e]0, +oo[xR 3 x]0, +oo[. 

We refer to p[5] for a global existence theorem for the discrete diffusive coagulation - 
fragmentation equation in L , and to |17] for the continuous one, improved in [19] (with 
less restrictive conditions on the kernels), then in [2] (with uniqueness of the solution). 



The second way to correct the spatially homogeneous problem is to assume that 
the particles are transported with a deterministic velocit y v. At the statistical level, 
this adds a linear transport term v.V x f to the equation (11.31) . This velocity can be a 
given velocity v = v(t,x,m) or the inner velocity of the particles. The first case has 
been studied in [6], with an existence and uniqueness theorem, and furthermore the 
continuous dependance on the initial data. Physically, it corresponds to the dynamics 
of particles with rather low mass which follow a velocity drift depending only on the 
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surrounding fluid. In the second case, particles are also identified by their momen- 
tum p G M 3 in addition to their mass m (with v = p/m): we have a kinetic model, 
which is relevant to describe the dynamics of particles of varying size/mass according 
to coagulation/fragmentation events, like in aerosols. At the microscopic scale, the co- 
agulation/fragmentation processes become "multi-dimensional", with mass-momentum 
conservation at each merger/break up according to the following scheme: 

Coagulation : {m} + {m*} — > {m'} 
{p} + {p*} -> {p>} 

Fragmentation : {m!} — > {m} + {m*} 

{p} - M + {p*} 
where m! := m + m*, m > 0, m* > 0, and p' := p + p*. 

Thus, in the statistical description, the density depends on time, position, mass and 
momentum: / = f(t,x,m,p). But even if this kind of kinetic models provides a rather 
good description of phenomena, it is harder to study, so there are less results than 
for the diffusive ones. Moreover, it is difficult to know the exact physical form of 
the kernels. And finally, the numerical aspects are a real problem on these models: 
because of a high dimension (at least 7 plus time), it seems to be very difficult, maybe 
impossible, to compute the solutions on a long time. 

Concerning the results, a global existence theorem for the sole coagulation has 
been demonstrated in [7]. The proof is based on L p -norms dissipation for any formal 
solution, and on weak-compactness methods in L . This result has been extended to 
a more general class of coalescence operators in [12j (but under stronger restriction on 
the initial data), with a very different method of proof. For the sole fragmentation, 
a difficulty is due to the blow-up of kinetic energy, which grows at each microscopic 
event. Thus, it is reasonable to take the internal energy of particles into account, which 
bala nces the gain of kinetic energy during a break up. With that modeling, the work 
provides global existence for a kinetic fragmentation model, with general growth 
assumptions on the kernel B, by using correct entropies. 

The aim of this work is to combine both of these analysis. The main difficulty is 
the lack of a priori estimates, because of the compensation problems between the two 
kernels: we are not able to define well the reaction term of coagulation with only the 
a priori bounds (specifically, we can not say that Qc(f,f) lies in L loc ). That is why 
we use the DiPerna-Lions theory of renormalized solutions, introduced in [3] to show 
global existence for Boltzmann equation, which presents the same problems. 



1.3 Description of the kinetic model and outline of the paper 

Now, let us describe precisely the model we study. The parameters which describe 
the state of a particle are denoted by 

y := (m,p,e) G Y :=]0, +oo[xIR 3 x]0, +oo[, 

m for the mass, p the impulsion, and e the internal energy. At the microscopic scale, 
coalescence and fragmentation conserve total energy (kinetic energy + internal energy), 
thus we can compute the internal energy of daughter(s) particle(s). 

Coagulation : {e} + {e*} — > {e} 
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We have 

g + e+ bl! + e *= \P + ^\ +e >, 
2m 2m* 2(m + m*) 

thus 

/ ^ / 4. 4.n , ^ / 4. |m*p — mp*\ 2 

e — e + e + E-(m,m ,p,p ), where E_(m,m,p,p):= — ^ > 

2mm* (m + m*) 

(E- is the loss of kinetic energy resulting from the merger). 

Fragmentation : {e'} — > {e} + {e*} 

We have 

\pf , \p\ 2 \p'-p\ 2 

— — he = — — h e H — — — h e 

2m' 2m 2(m' — m) 

thus 

e = e — e — E + (m ,m,p ,p), where EAm ,m,p ,p) := -, — > U 

2mm\m — m) 

(E + is the gain of kinetic energy resulting from the break up) . 
Remark 1.1 Let us point out the following symmetries: 

E-(m*,m,p*,p) = E-(m, m*,p, p*) and E + {m ,m,p ,p) — E + (m! ,m! —m,p' ,p' —p), 

and the relation: E_(m,m* ,p,p*) = E + {m + m* , m, p + p* , p) , which is consistent with 
the two phenomena's reciprocity. 



We use the following notations: 

• if y = (m,p,e), y* = (m*,p*,e*), then we denote 

y := y + y* := (m + m* ,p + p*, e + e* + E_(m, m*,p,p*)), 

• if y = (m,p,e), y = (m',p',e), with m < m and e < e — E + (m ,m,p ,p), then 
we say that y < y' and we denote 

V* -=y' ~y ■= (m' - m,p' - p,e' - e - E + (m' \m,p' \p)). 

With this formalism, we naturally have {y —y)+y — y , but note carefully that y <y 
is not an order relation on Y . 

Remark 1.2 For ally 1 E Y, {yeY, y< y'} c]0, m'[xB ^^^ x}0, e'[. 
Denoting Y R :=]0, R[xB R x]0, R[c Y, we have 

y<y', y'EY R y e]0, R[xB^ R x]0, R[c Y m . (1.5) 
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Finally, we point out that the map {m! , m, p' , p, e' , e) i— > (m', m*,p',p*, e', e*) is a dif- 
feomorphism with C°° regularity whithin the domain 

{0 < ra < to', G M 3 , < e < e' - E + (m',m,p',p)} C F 2 
which preserves volume. 

We denote by f(t,x,m,p,e) = f(t,x,y) the particles density, which is a nonnega- 
tive function depending on time t > 0, position x G M 3 , and the mass-momentum- 
energy variable y. To shorten the notations, we set for each t,x, f(y) = f(t,x,y), or 
/ = f(t,x,y), f* = f(t,x,y*), and /' = f{t,x,y ). The complete model then reads: 



p 



d t f + ^.V x / = Q c + (/, /) - g c -(/, /) + Qj(/) - QJ(/), (£CF) 



t E]0, +oo[, iGl 3 , y = (m,p,e)eY, 



with 



and 



' Qc (/> /)(v) = 5 jf ^(v*, y - v*)f(y*)f(v - y*)^< y} dy\ 

Qc(fJ)(v) = f(y)Lf(y), Lf(y) := / A(y,y*)f(y*)dy*, 
<#(/)(v) = J Y B{y',y)f{y')t w>y} dy', 

Qj(f)(y) = lB 1 (y)f(y), B 1 (y') := J B(y', y)l {y<y , } dy. 

Functions A et B are respectively the coagulation and fragmentation kernels. They 
are nonnegative functions, independent of (t,x), which satisfy the natural properties 
of symmetry: 



V(y,y*)eY 2 , A(y,y*)=A(y\y), 



(1.6) 



V(y,y')eY 2 , y<y', B(y',y) = B(y',y*). (1.7) 

The kernel A(y,y*) represents the coalescence rate between two particles y and y*, 
whereas B(y' , y) is the fragmentation rate for a particle y' which breaks in two clusters 
y and y*. 



We assume that A fulfils the following structure assumption: 

V(y, y*) G F 2 , A(y, y*) < A(y, y') + A(y*, y')- 
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Remark 1.3 We can insist on the fact that this assumption is more general than the 
classical Galkin-Tupchiev monotonicity condition: 



Vy<y\ A(y,y*-y)<A(y,y*). 



;i.g) 



In th e "monodime nsion al" case, the Smoluchowski kernel given by do not satisfy 
d 1 . 9i) but satisfies f 1 1 . <S1) . that's why the first existence result established in flTtf under 
Galkin-Tupchiev condition was extended in [WJ to kernels which satisfy f 1 1 . 8i) only. 



We also require that A and B have a mild growth: 



VR > 0, 



Vi? > 0, 
and B is truncated as: 



y r \W 
B(y',y) 



dy — > 0, 

|y*|^+oo 



Yr 



l {y<y > } dy — » 0, 
y | |i/'|-H-oo 



(1.10) 
(1.11) 



3C n > 1, 



or 



m! > C m 



2m' ° \ ^ 2m 



B(y\y)=0. 



;i.!2) 



Remark 1.4 The physical interpretation of this truncature assumption is to prevent 
the creation of too small clusters compared to the mother particle. From a mathematical 
point of view, it allows the total number of particles (the L l -norm of f) to be finite at 
each time t > 0. 



We also need to have B\ locally bounded: 

V J R>0, B.eL^iYn) 



as well as A: 



V J R>0, A G L°°(Y^). 



;i.i3) 
;i.i4) 



Remark 1.5 Unfortunately, these assumptions of growth and boundedn ess a re more 
restrictive, and in the monodimensional case, the Smoluchowski kernel (Q3) doesn't 
satisfy them any more. The examples given in [2\j for the sole coagulation, namely 



A(m, m*, p, p*) = (m a + m* a ) 2 



p p 



0<a< 1/2, 



(for the dynamics of liquid droplets carried by a gaseous phase) or 



A(m, m* ,p,p* 



fm + m* N 


f 


P_ _ 




\ mm* / 




m 


m* 



< a < 1, -3 < 7 < 0, 



(for a stellar dynamics context) do not fit neither. Here, we need coalescence kernels 
which are bounded when m,m* — > 0. But it is difficult to know the exact physical form 
of the kernels A and B because of the complexity of this kinetic model. Nevertheless, 
simple kernels given by A(m, m*) = m a + m* a with < a < 1 fit. 
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Finally, we assume that A controls B in the following sense: 

3s > 1, 30 < 5 < — - — < 1, 
6s — 5 

Remark 1.6 This last assumption is more technical, but seems necessary to balance 
the contributions of the interaction terms Q c (f,f) and Qf(f), which are difficult to 
compare because Q c (f,f) is quadratic whereas Qf(f) is linear. 

The paper consists in the proof of the following theorem. 

Theorem 1.1 Let A and B be kernels satisfying ( U.bA) — and (EZ/j) — (QU3) and 

let f be a nonnegative initial data which satisfies 

1 + m + — — h e + m\x\ 2 J f°(x, y) + f°(x, y) s ) dxdy < oo, 
2m / / 

(1.16) 

then for all T > 0, there exists f E C([0, T], L X (R 3 x Y)) such that f(0) = f° and f is 
a renormalized solution to (ECF). Moreover, 

a.ete]0,T[, [[ (l+m+^ + e + mlxl 2 ) f(t,x,y)dxdy<K T , (1.17) 

o.ete]0,r[, // f(t,x,y) s dxdy < K T , (1.18) 
where the constant Kt depends only on Co, T, K(f°), s and S (defined in <\1. and 

Beyond existence problems, there are lots of others interesting subjects to explore. A 
first one concerns the mass conservation of the solution /, which is still an open problem 
for such kinetic models, even for the case of the sole coagulation. In the spatially 
homogeneous case, it has been shown in [5] that total mass is preserved in time under 
mild growth hypotheses on kernels. But we know that in case of strong coagulation 
(typically the case of multiplicative kernels), a phenomenon of gelation occurs, which 
force the total mass of the system to decay from a certain time T g < +oo. Then, 
problems of convergence to an equilibrium have been already studied for the spatially 
homogeneous equation [H], under a detailed balance condition between kernels A and 
B. We can also mention existence of self-similar solutions [H M, E], always for the 
spatially homogeneous case. 

In a first section, we will derive the a priori estimates from the equation, giving the 
proper setting of the problem. Then, the proof of theorem is based on a well-known 
stability principle which says that if we are able to pass to the limit in the equation (the 
set of solutions is closed in a certain sense), then it would be easy to show the existence 
of a solution, applying the stability result to a sequence of approached problems which 
we can solve. So, the aim of the last section is to prove rigorously such a stability result 
and in fact, we work in the context of renormalized solutions, because the reaction term 
can not be defined as a distribution simply using the a priori estimates. 
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1.4 Different notions of solutions 



We discuss here on different notions of solutions, recalling the DiPerna-Lions results. 
We set Q(f, f) = Q+(/, /) - Q~{f, f) + Q+(f) - Qj(f). 

Definition 1.2 Let f be a nonnegative function, such that f G L^QO, +oo[xIR 3 x Y). 
We say that f is a renormalized solution of (ECF) if 

O^fl G L l oc(]0) +oo[xR 3 x r)? QfSV G LL(]0) +oo[xR 3 x n 

and if the function g := log(l + /) satisfies the renormalized equation 

d t9 + ^-V x g = (ECFR) 
m 1 + / 

in D'CO^oofxR 3 x Y). 

The renormalization makes passing to the limit impossible because of the quotients 
in the reaction term, that is why we also need another notion of solution: the mild 
solutions, which only require local integrability in time and provide Duhamel's integral 
formulations to the problem in which we are able to pass to the limit. 

Definition 1.3 Let f be a nonnegative function, such that f E L loc (]0, +oo[xIR 3 x Y). 
We say that f is a mild solution of (ECF) if for almost all (x,y) G M 3 xY, 



VT > 0, Qt(f,f)Kt,x,y) G L\]0,T\), Qf(ff(t,x,y) G L l (]0,T[), 

and 

V0 < s < t < oo, f% x, y) - f(s, x,y) = J Q(f, ff(a, x, y)da, (1.19) 

where hs denotes the restriction to the characteristic lines of the equation: 

hfi(t, x, m,p, e) := h(t, x + t — , m,p, e). 

m 

The following results are proved in [3j: 
Lemma 1.4 

(i) IfQt(fJ) e LL(]0,+oo[xM 3 x Y) and Qf(f) G LL(]0,+oo[xM 3 x Y), then 
the following assertions are equivalent: 

• f is a solution of (ECF) in the sense of distributions. 

• f is a renormalized solution of (ECF). 

• f is a mild solution of (ECF). 
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(ii) If f is a renormalized solution of (ECF), then for all function (3 G C 1 ([0, +00 [) 
C 

such that \f3'(u)\ < the composed function (3(f) is a solution of 

l + u 

d t (3(f) + ?-.V x (3(f)=P'(f)Q(fJ). 
m 

in the sense of distributions (here, the right side lies in Lj oc (]0, +oo[xM 3 x Y) ). 

(Hi) f is a renormalized solution of (ECF) if and only if 

Q±( f f) Q ± (f) 
f is a mild solution of (ECF), [ ^ e L L(1°> +°o[xM 3 x Y) and y^J E 

Ll c (}0,+oo[xR 3 xY). 

2 A priori estimates 

We consider the Cauchy problem 

{ECF) (2 20) 

f(0,x,y) = f°(x,y). {Z - ZU) 

We suppose in this section that (I2.20p admit a sufficiently smooth solution / in order 
to handle some formal quantities which are conserved or propagated by the equation 
(ECF). More precisely, we will show the propagation of L q bounds for the solution 
along time: 

Proposition 2.1 If the initial data f° satisfies 



bl 



2 



K(f°) :=// 1 + m + ^ + e + m\x\ 2 f°(x, y) + f°(x, y) s dxdy < 00, 



2m 

(2.21) 

then for all T > 0, any classical solution of the Cauchy problem \2.20) satisfies 



sup // ( ( 1 + m + ^ + e + mix] 2 ] f(t,x,y) + f(t,x,y) q ] dxdy < K 
m,T]JJm?xY \\ 2m J J 



te[o,T] 

for all the exponents q g] 5/6, s], and also 



T- 

(2.22) 



Ft (D 1 (f(t,x)) + D 2 (f(t,x)))dxdt<K T , (2.23) 
Jo Jr 3 



where 



D 1 (f(t,x)):=l II A(y,y*)sup(f,r)M(f,r) s dy*dy>0, (2.24) 



D 2 (f(t,x)) : = i-i / / B(y',y)f(t,x,y') s % <y r } dy'dy > 0, (2.25) 
and the constant depends only on Co, T , K(f°), s and 5. 
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2.1 Basic physical estimates 



We start with a fundamental formula, which gives the variation in time of some integral 
quantities involving the solution /. 

Lemma 2.2 Let H{u) be a function with C 1 regularity on [0, +oo[ and $(y) a real or 
vectorial function. We have 



d 
di 

1 

2 



If $(y)H(f(t,x,y))dxdy 

Aff* (&d u H(f) - $d u H(f) - &d u H(t)) dy*dydx 




<xYxY 



+ l Iff B f (® d u H (f) + ®*d u H(f*) - <f>'d u H(f')) U {y<y , } dydy'dx, 

1 J J JR^xYxY 

(2.26) 

where d u H = — — . 

du 

Proof : Using (ECF), we have 

4// <f>(y)H(f)dxdy = ff $(y)d u H(f)d t fdydx 

= ff S(y) d u H(f) (Qt(f, f) - Q:(f, /)) dydx 

J JR 3 xY 

+ ff $(y)d u H(f) (Q+U)-Q- f U))dydx 

J JR 3 xY 

- ff div x (-<S>(y)H(f)?-)dydx. 

J JR^xY m 

The integral with divergence vanishes thanks to Stokes' formula. Whence 

± [f ®{y)H{f)dxdy = ff $(y) d u H(f) (Q+(f, f) - Q~(f, /)) dydx 

Q>T J JR^xY J JR z xY 

+ 11 $(y)d u H(f) (Q+(f)-Q](f))dydx. 

J JR 3 xY 

Using Fubini's theorem (formally), we can write 
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d 
~dt 





'■xYxY 



•xYxY 



+ 



[ [ $(y)H(f)dxdy 

J JR 3 xY 

$(y) d u H(f) A(y\ y - y*)f(y*)f(y - y*)l {r<y} dy*dydx 
$(y) d u H(f)A(y, y*)f(y)f(y*)dy*dydx 
[[[ $(y) d u H(f) B(y',y)f(y')l w>y} dy'dydx 

J J JRZxYxY 

-\ff ^d^U^B^fWy'dx. 

z J JR 3 xY 

If we change variables (y* ', y — y*) — > (y*, y) in the first integral, we obtain 
± [[ <S>(y)H(f)dxdy 

al J JR 3 xY 

$(y + y*) d„ff (/(y + y*)) A(y*, y)f(y")f(y)dy^dydx 
$(y) d u H(f)A(y, y*)f(y)f(y*)dy*dydx 




'•xYxY 




'•xYxY 




+ /// $(y)^(/) J B(y',y)/(y / )%'>,}^'^ 




$(y') d tt if (/') 5(y', y)/(yO% <2/}( W^. 

; x7xy 



The symmetry of A allows us to write 

$(y) d u H(f)A(y, y*)f(y)f(y*)dy*dydx 




•xYxY 




$(y) d u H(f)A(y, y*)f(y)f(y*)dy*dydx 

h xYxY 



xYxY 




2 

using the change of variables (y, y*) — > (y*,y). 
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The same applies to B with (y',y) — * (y 1 , y' — y): 

fff <S>(y)d u H(f)B(y\y)f(y')l {y , >yy dy'dydx 

J J JR 3 xYxY 

<S>(y)d u H(f)B(y',y)f(y')l {y , >y} dy'dydx 




'•xYxY 



+ l I I I Hy' -y)d u H(f(y' -y))B(y',y)f(y')l {y , >y} dy'dydx. 

£ .1 .1 ./w3 x y x y 




□ 



Applying this lemma with H(u) = u, it gives 

11 Hy)fdxdy = I HI Aff ($'-$-$*) dy*dydx 



dt J ./d3vV 2 



xY * J J JR?xYxY 



+ l 1 1 1 Bf($ + &- $') l {y<y/} dy'dydx. 

"'XYxY 




(2.27) 



2 

Choosing $(y) = to, we obtain mass conservation: 

^- ff mf(t,x,y)dxdy = 0. (2.28) 

at J JR 3 xY 

With $(y) = p, we also get the momentum conservation: 

ff pf(t,x,y)dxdy = 0. (2.29) 

J JR 3 xY 



d_ 

~dt 



\v\ 2 

Then, choosing $(y) = h e, we recover the total energy conservation: 

2to 

jtILM +e ) nt - x ' y)dxdy=o - (2 - 30) 

Moreover, we can control space momenta: 

Lemma 2.3 For all T > 0, there exists a constant Ct > that 

Vte[0,T], // m\x\ 2 f(t,x,y)dxdy <C T . (2.31) 

J JR 3 xY 
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Proof : In view of the equation (ECF) and the Stokes formula, we have 



d 
di 



m\x\ 2 fdxdy 



<xY 



\x\ 2 p.V x fdxdy 



<xY 



x.p f(t, x, y)dxdy 

<xY 

v 1/2 

m\x\ 2 fdxdy J 
and we conclude with (12.301 ) and Gronwall's lemma. 



< 2 



\JJr3xy m J 



1/2 



□ 



Finally, we can control the number of particles in finite time: 
Lemma 2.4 We set 



N :-- 



En :- 



f°(x,y)dxdy, M :-- 



•xY 



'•xY 



mf°(x,y)dxdy, 



■xY 



+ e ) f°(x. !j)(l.nhj. 



Then, there exists a constant C > depending only on Cq that 

VT > 0, Vi G [0, T], II f{t, x, y)dxdy < (N + CT(M + E ))e CT + M + E . 

J JR 3 xY 

(2.32) 



Proof : We use formula (12.271) with $(?/) = 1, <t g b£ <1 i- Since $ is nonnegative and 
subadditive in the sense of coalescence (ie $' < $ + $*), we have 



d 
d£ 



1, ^ i„|2 ^ fdydx < - 

{™<1, e+fL<l}' / y - 2 




! xyxr 



5/' ($ + $*- $') ]l {y<y , } d ? /d|/ / da; 




! xrxY 



S/' ( $ - y ) t {y<y r } dydy'dx 




< I I I Bf'$l {y<y , } dydy'dx 

'xYxY 




•xYxY 



Bf'l. , „ |„|2 ^dydy'dx. 
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In the last integral, if m' > Co, then, since m < 1, we have B(y',y) = according to 

\P'\ 2 

assumption (11.121) . The same applies if e H > Co. Thus, 

2m' 

-r ff f(t,x,y)t r „ i„|2 ^dydx 



< 



'xy \jy 



B(y',y)l {y<y , } dy j f{t,x,y')l {m ,^ Coel+ ^ Co} dy'dx 



Denoting C := sup Bi(y'), we obtain 
~r ff f(t,x,y)l c ^ |„|2 ^dydx 

dtJJjexY {m<l,e+M.<l} y 



< c 1 1 f{t,x,y)dydx 

■xY 



< c 



f(t,x,y)l | P |2 dydx+ // mf(t,x,y)dydx 

II ( e +7r~) f(t,x,y)dydx] . 



'xY 



2m J 



Using (I2T28D and ( 12301 ). we have 

R 3 xY 



/(*^^)V<l,e+r<l}^ 

L ' 2 171 — J 



- c yy R3xy ^ x ^)v<i,e + ^<i}^ x+c ( Mo+jE;o )- 

We integrate this inequality in time. Then, Gronwall's lemma provides 

VT>0, Vte[0,T], // f(t,x,y)t M * dydx<(N +CT(M +E ))e 

We conclude noting that 

// f(t,x,y)dydx < // f(t,x,y)l M 2 dydx 



^m + e+ f(t,x,y)dydx, 
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and using (12 . 281) and (12.301) again. 

□ 

\p\ 2 

To summarize, if we set E(x, y) = 1 + m H h e + mlrd 2 , and if we suppose that the 

2m 

initial data satisfies 

K{f) := !! E(x, y)f°(x, y)dxdy < +oc, 

J JR 3 xY 

then, for all T > 0, there exists a constant i^r (depending on T, Co and K(f )) such 
that 

sup // E(x 1 y)f{t 1 x,y)dxdy<KT- (2.33) 

te[0,T]J Jr^xY 

Remark 2.1 For 7 > 5, we have 

1 



M 3 xy EF(x,y) 



dxdy < +00. (2.34) 



It tui/Z 6e ven/ useful to show that some L q bounds of f (for 5/6 < q < 1 and q = s > 1 ) 
also propagate in time. 



2.2 L q bounds 

Obtaining L q bounds propagation is more technical, that is why we split the proof in 
several lemmas. 

Lemma 2.5 Let (3 6] 5/6, 1 [. Then, for all T > 0, there exists a constant (de- 
pending on T , Co and K(f°) ) such that 

sup // f(t,x,y)dxdy < K T . (2.35) 

t&[0,T] J JR^xY 

Proof : Writing 

f /3 (t,x,y)dxdy= [[ fP(t,x,y) E ^' V \ dxdy, 
■xy JJrSxy EP(x,y) 

we use Young inequality 

u a v a * 

Vet > 1, Vu > 0, W > 0, uv < — + (2.36) 

a a* 

with u = f (3 (t,x,y)E l3 (x,y), v = — ^ — r, and a = - > 1, 

EP(x } y) (3 
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and obtain 

f /3 (t,x,y)dxdy < (3 E(x,y)f(t,x,y)dxdy 

<xY J JR 3 xY 

+ (1-/3) If dxdy. 

J Jr3 xY E*-p (x, y) 

We conclude with JZHSU and (I234D . 

□ 

Lemma 2.6 For any convex and nonnegative function H E C 1 ([0, +oo[) such that 
H(0) = 0, and for all t > 0, we have 



[[ H(f(t,x,y))dxdy < [[ H(f°(x,y))dxdy 

J JR s xY J Jr 3 xY 



fllf A sup(/, f)H(mf(f, F)WdydxdT 

loJ J JR^xYxY 

A'H f'A {y<yl} dydy'dxdT 
Bfd u H(f')t {y<yl} dydy'dxdr, 




OJ J JR :i xYxY 



1 '* 




+ 



2 JQJ J jR^xYxY 

(2.37) 

where A = A{y,y*), A' = A(y,y'), B = B{y',y). 
Proof : The formula (12.261) yields 

^-11 H{f)dxdy = Iff Aff (^Eill _ d u H(f)) dy*dydx 

at J JR?xY J J JRZxYxY \ Z J 



+ [[[ Bf (d u H(f) - ^Eill) t {y<yl} dydy'dx. 

J J JR?xYxY V 1 J 



(2.38) 

Let us rewrite the term I\ := / / / Af f*d u H(f')dy*dydx, by the following way: 

^xYxY 




h= A inf (/, f) sup(/, nd u H(f')dy*dydx. 

J J JR^xYxY 

We use the Young inequality: 

Vu > 0, > 0, < H(u) + (2.39) 



17 



with u = sup(/, /*) and v = d u H(f), where H* stands for the convex conjugate 
function of H. A simple calculus shows that 

H*(d u H(u)) = ud u H{u) - H(u), 

and this quantity is nonnegative, by the assumptions on H. We denote 

0(u) := H*(d u H(u)) > 0. 

It leads to the inequality: 

h < III Aini(f,f*)H(suv(f,r))dy*dydx 

J J JR 3 xYxY 





+ /// AM(f,r)Q(f)dy*dydx. 

J J JR 3 xYxY 

We can dominate the second term of the right member using the hypothesis (11.81) by 
the following way: 



Amf(f,r)G(f')dy*dydx 

•xYxY 

< I 1 1 (A(y, y + y*) + A(y\ y + y*)) inf (/, fW)dy k dydx 

'■xYxY 

< HI A(y, y + y*)fQ(f')dy*dydx + /// A{y\ y + y*)f*e(f')dy*dydx 

'xYxY J J JR 3 xYxY 

2 / if A{y\ y + y*)/(y*)6(/(y + y*))dy*dydx 

! xYxY 

A(y\y)f(y*)e(f(y))t {y * <y} dy*dydx 

s xYxY 

(the last identity resulting from the change of variables (y*,y + y*) — > (y*,y)). 
This yields 

h < HI AM(fJ*)H(suv(fJ*))dy*dydx 

J J JR 3 xYxY 






+ 2 Af*e(f)l {y * <y} dy*dydx. 

J J JR^xYxY 

Thus, we have the following control of the coagulation contribution in (|2.38l) : 
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J J Jm 3 



! xYxy 



< \ I I l Ainf(/, (sup(/, r))dy*dydx 

'•xYxY 





■xYxY 



2 

+ /// Af*Q(f)l {r<y} dy*dydx - III Afrd u H(f)dy*dydx. 

J J JR 3 xYxY J J JR 3 xYxY 

Now, we can write the right member of this inequality by the following way: 

Asup(fJ*)H(mi(fJ*))dy*dydx 
I [[[ AfH(f*)dy*dydx + \ [[ [ AfH(f)dy*dydx 

1 J J JWxYxY 1 J J JW'xYxY 

- / / / Af*e(f)l {r<y} dy*dydx - [[[ Aff*d u H(f)dy*dydx 

xYxY J J JR^xYxY 

l [[[ Asup(/, (inf(/, t)Wdydx 

Iff Ar<d(f)l {r<y} dy*dydx -Iff Are(f)dy*dyd 

J J JR^xYxY J J JR s xYxY 





+ III \rU(t)l U r. „\<l!/ <'!/<!■>■ Ill Aj^-)(j)(l„(l,)(l,V. 

We deduce 



/// Aff (^P- - d u H(f)) dy*dydx 

J J JR 3 xYxY \ 1 J 

< \\\\ ^sup(/, f)H(ini(f, nWdydx 

1 J J JR^xYxY 



+ /// AfQ(f)l {r<y} dy*dydx - [[[ AfQtfWdydx. 

J J Jr^xYxY J J JR s xYxY 

(2.40) 

Then, we can also control the fragmentation contribution: 

Bfd u H{f)l {y<y , } dydy'dx -Iff Bf^p-l {y<y , } dydy'dx. 

l xYxY J J Jr^xYxY 1 

We rewrite the first term by the following way: 
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Bfd u H{f)t {y<y , } dydy'dx = \\\ ^A' f'd u H(f)l {y<y , } dydy'dx 



"■xYxY J J JR A xYxY 



and we use (12.391) again, with u = — and v = d u H(f), whence 




■xYxY 



< 




Bf" fd u H(f) - duH (f ^ ^ l {y<y/} dydy'dx 

H ( — j A'f'l {y<y , } dydy'dx + [ [ [ A' f'Q(f)l {y<y , } dydy'dx 

\ A J J J JR 3 xYxY 



3 xYxY 




Bf' duH ^ t{ y<y > } dydy'dx. 

■xYxY 1 

(2.41) 

Eventually, using (ESH)), li2T4"Gl and (EHj), we infer 

^ // H(f(t,x,y))dxdy < \lll Asup(/, t)H(mf(f, r))dy*dydx 

at J JR 3 xY A J J JR 3 xYxY 

Af*Q(f)dydy*dx 

'■x{Y2-({y<y*}U{y*<y})) 

+ Iff A'H f'l {y<yl} dydy'dx 

J J JR 3 xYxY V^ 1 / 

~~ III B f' duH y - \ y< y>}dydy'dx. 

J J JR 3 xYxY 1 

(2.42) 
□ 

Lemma 2.7 Forall T > 0, there exists a constant Ct > depending only on T, the 
initial values N ,M ,E and the truncature parameter Cq such that for all t G [0, T], 




// f(t,x,y) s dxdy < // f°(x } y) s dxdy + C 7 

J JR 3 xY J JR 3 xY 



nil A sup(/, f*) inf (/, n s dy*dydxdr 

IQ J J JR 3 xYxY 



s-S 



B (f') s l {y<yl} dy'dydxdr, 

'xYxY 




(2.43) 

where s and 5 are given by $1.150 . 
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Proof : We use the previous lemma with H(u) = u s . We obtain 

f(t,x,y) s dxdy < II f°(x,y) s dxdy 



I "' 
2 




A sup(/, /*) inf (/, f*) s dy*dydxdr 

oJ J Jr 3 xYxY 



' B 




J il 3 xYxY 




f'l{ y< y'}dydy'dxdr 
B 1 (y')(f , ) s dy'dxdT. 



According to (fTTTBl) and fl2~33l) . 




< 



) f'\ y <y'}dydy'dxdT 




JoJ JR 3 xY 

1 " 



(\ +rri + ^ + «A fdy'dxdr + - f ff B 1 (y , ) s fdy'dxdr 
V 2m ' / ZJoJJrsxy 



< TK T + - B 1 ( y ') s f' sS f' 1 - sS dy'dxdr. 
1 JoJ Jr 3 xY 

We apply the Young inequality again with the exponent 1/5 > 1: 

{B,{y'Yr s ) ll& (/"-*) <W* 



1/5 (1/*)* 

1 — s<5 

5 J B 1 ( 2 / / )r + (i-5)/'— . 



Thus we deduce 



// f(t,x,y) s dxdy < // f°(x,y) s dxdy 
J Jm 3 xY J Jr 3 xy 



2 




A sup(/, /*) inf (/, f*) s dy*dydxdr 

IqJ J Jr 3 xYxY 



+ TK T + iy- / / / f'^dy'dxdr 




OJ JR 3 xY 



fff B 1 (y')(fydy'dxdT. 
1 JoJ Jr 3 xy 
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. . 1 — s8 , 

We can use (E2D) since £ 5/6, 1 . 

l — o 



f(t,x,y) s dxdy < // f(x,y) s dxdy + C T 

<xY J JR 3 xY 

1 " 



2 




Asup(/, /*) inf(/, f*) s dy*dydxdT 

OJ J JRfixYxY 



fit B^if'Ydy'dxdr. 



We conclude noting that 5 < < 1 < s. 

6s — 5 



□ 



3 A stability result 

The proof of theorem [TTT1 relies on a stability theorem, which claims that we can pass to 
the limit in the equation (ECF) in a certain sense, namely in an integral formulation. 

Definition 3.1 Let T > and let f be a nonnegative initial data which satisfies 
\2.21\) . A weak solution of [2.2(A) is a nonnegative function f G C([0, T], L 1 (M 3 x Y)), 
verifying the estimates $2.22]) and $2.23\) , satisfying (ECF) in £>'(]0, +oo[xR 3 x Y), 
and such that f(0) = / . 

Now let us state the result we will prove in this section: 

Theorem 3.2 Let (f n )n>i be a sequence of weak solutions of [2.2(h) . with initial data 
/°, and such that 

\/n G N, fn e W l ' l (]0, +oo[xM 3 x Y), (3.44) 
sup sup ff ( ( 1 + m + ^— + e + m\x\ 2 J f n (t,x,y) + f n (t,x,y) q J dxdy < K T , 

n>l te[0,T] J JR3 xY V V -* m J J 

(3.45) 

for all the exponents q G] 5/6, s], and also 

sup / / (D 1 (f n (t 1 x)) + D 2 (f n (t, x))) dxdt < K T . (3.46) 

n>l Jo JR 3 

(the a priori estimates hold uniformly in n). 

Then, up to a subsequence, f n —^f weakly in L 1 (]0, T[x]R 3 oc x Y), where f is a renor- 
malized solution of (ECF). Furthermore, f G C([0, T],L 1 (R 3 x Y)). 
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3.1 Weak compactness of (f n ) 

Let < T < oo. The bounds on f n provides some weak compactness, and thus the 
existence of a limit / after extraction. 

Lemma 3.3 For all R > 0, the sequence (f n )n>i is weakly compact in -^ 1 (]0, T[xB R x 
Y). 

\p\ 2 

Proof : We set := f s and \& (m.p.e) := m H h e. The function $ is nonde- 

J v ' 2m 

creasing, nonnegative and $(£)/£ — > +°°) ^ is nonnegative and ^(y) — ► +oo. 
The estimate (13.451) gives 

sup/ / / ((l + tf(y))/ n + $(/ n ))dtdxdy<+oo, 



n>l JO JBrJY 

and we conclude by Dunford-Pettis theorem. 

□ 

Thus, there exists a nonnegative function / such that for all R > 0, f n — / in 
L 1 (]0, Tfxi?^ x F) for a subsequence (not relabeled). Moreover, we can show easily 
(diagonal extraction) that the subsequence is not depending on R. Then we notice 
that in fact, / G L°°(]0, T[, L\R 3 x Y)) and 



\P\ 
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a.ete]0,T[, // 1 + m + ^- + e + m\x\ 2 ) f(t,x,y)dxdy < K T . (3.47) 
Moreover, since the function £ i— > is convex, we have 

a.ete]0,T[, // f(t,x,y) s dxdy < K T . (3.48) 

3.2 Weak compactness of the renormalized coalescence term 

The bounds on f n are not enough to define the term Qc(fmfn) as a distribution, 
unlike the term Qt(fmfn)- So, it seems that renormalization is necessary to obtain 
well-defined and weakly compact coalescence terms. 

Lemma 3.4 For all R > 0, the sequence (Qtifn, /n))n>i is weakly compact in 
L x (]0,T[xB R x Y R ), where Y r :=]0,R[xB r x]0,R[. 

Proof : Let E be a measurable subset of ]0, T[xB R x Y R . We set ip(t, x, y) := fl^i, x, y). 
Performing the change of variables (y,y*) — > (y*, y — y*), we obtain 
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QUfn, fn){y)^{y)dy = \ I I A(y*, y - y*) f n (y*) f n (y - y*)(p(y)l {y * <y} dy*dy 



'Yr 2 JYrJY 





1 

2 

0<m<R 0<m*<R-m 

peM 3 p*eB(-p,R) 

0<e<R 0<e*<R-e-E-(m,rn* ,p,p*) 



In fact, we only integrate over p G B 2 r because 



My, y*)fn(y)fn(y*My + y*)d y *d y . 



2m ~ 2m 2m* 2(m + m*) v ' '^'^ ; ~ 2(m + m*) 



which yields 



U|2 < ^2^^ + 2mjR < 3jR 2 

m + m* 



(and the same applies to p* because of the symmetry in the previous computation). 
Thus we have 



f Y Qt(fnJ n ){y)v{y)dy = \ JJJ JJJ Af n ftfdy*dy. 



0<m<R 0<m*<R-m 
peB 2R p*€B(~p,R) 
0<e<R 0<e*<R-e-E-(rn,m*,p,p*) 



Using the inequality 

Af n f* < j^A S Up(f n , /*) inf(/ n , ^) S %„f( /n ,/*)>M} + MAsU P (/„, /*)fl{inf(/„,/*)<M}, 

(3.49) 

we obtain 

f QUU, fn)(yMy)dy < + T / / ^sup(/„J n >'^y 

Jy r M s 2 J Y2R Jy 2R 



DMt.x)) 
M s 



Y 2 r JY 2 r 



< 



D!(f n (t,x)) 



M* 
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We fix e > and choose M such that < e. 

So, we can write 



Qt(fnJn)(yMy)dy < sD^f^x)) + M^IU^ ■ / f n v'dy*dy 
Yr Jy 2R Jy 2R 



< eD x (f n {t, x)) + M\\A\\^ YiR / / f n <p*dy*dy. 




Y JY 



Eventually, in view of (13.451) and (13.461) . we obtain 

n\ Qt(fn,fn)tp{t,x,y)dydxdt < M||A|| 00i y 2 : / / / / f n tp*dy*dydxdt 
j R Jy r Jo Jb r JyJy 

+ eK T . 

We conclude by letting mes(E) — > and using the weak compactness of (/„). 



□ 

Corollary 3.5 For all R > 0, the sequence ( ^ ^ n, J n ^ ) ^ weakly compact in 

n>l 



l + /n 

^QO^lxfi^ x y fl ), w/iere Yr :=]0, i2[xS fl x]0, 



Proof : It's obvious by the previous lemma and Dunford-Pettis theorem since 

Qc jfnj fn) , /n+z' f f \ 
-i , r — \Jnj JnJ- 
J- ~r Jn 

□ 

Lemma 3.6 For all R > 0, Lf n — L/ weakly in L Q0,T[xBr x Yr), where Yr : = 
)0,R[xB R x)0,R[. 

Proof : Let <p(t,x,y) G L°°(]0, T[x5r x Yr). We have 



Lf n (yMy)dy= / / A(y,y*)f n (y*)cp(y)dy*dy. 

Y R JYrJY 

We fix £ > and, in view of the assumptions (11.61) and (11.101) . we choose R* > such 
that 

-4(2/, y* 



We can write 



Lf n ipdydxdt =1111 A(y,y*)f n (y*)tp(y)dy*dydxdt 
o Jb r jy r Jo Jb r Jy r Jy r * 



+ I [ [ [ A(y,y*)f n (y*) V (y)dy*dydxdt. 
Jo Jb r Jy r Jy-y r * 
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First. 



J Br JYrJYr* 



Af n (y*)tp(y)dy*dydxdt 



J B R JYrJYr 



Af{y*)v(y)dy*dydxdt. 



Indeed, setting 9(t,x,y*) = / A(y , y*)<p(t, x, y)dy , we have 



[ [ [ Af n {y*)^{y)dy*dydxdt= Iff 6(t,x,y*)f n (t,x,y*)dy*dxdt 
o ./Br JyrJYr* Jo J Br Jy r * 



and we conclude by lemma [3751 because the assumption (I1.14K implies 

0er(]OJ[xBiiX^). 
Moreover, 



J Br JYrJY-Yj 



Afn {y*)v {y)dy*dydxdt 



<£|MU/ / / \y*\f n {y*)dy*dxdt 
'o JbrJy-Yr* 



< eT H^IIoo-^t, 



and the inequality (13.471) yields 



Af(y*) V (y)dy*dydxdt 



io JBrJYrJy-Yr* 
Finally, we infer 



/ / Lf n tpdydxdt — I / / Lf (pdydxdt 
o J B R Jy r Jo J Br Jy r 



< o(l) +C(T,R,<p)e. 



□ 



Corollary 3.7 For all R > 0, the sequence 



Qc ifni fn, 
1 + fn 



is weakly compact in 



n>l 



L\]Q,T[xB R x Y R ), where Y R :=]0, R[xB R x]0, R[. 



Proof : It's obvious because ^ c (fn,fn) _ fn ^ 

1 ~r Jn 1 ~r /n 



□ 



3.3 Weak convergence of the fragmentation term 

Since they are linear, the fragmentation terms easily pass to the limit, and we have the 
following lemma. 

Lemma 3.8 For all R > 0, we have 
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(i) Qj(fn) ~± Q + f {f) weakly in L\]0,T[xB R x Y R ), 

(») Q](fn) Qf(f) weakly in L\}0,T[xB R x Y r ), 
where Y R :=}0, R[xB R x]0, R[. 

Proof : T he p art (ii) results immediately from (11 . 13H , and the proof of (i) is the same 
as lemma IcTBl 

□ 

3.4 Strong compactness of ^/-averages 

Strong compactness is needed to pass to the limit in coalescence terms (because they 
are quadratic), that's why we use the following averaging lemma, inspired by [7], [3], 
and [4]: 

Theorem 3.9 Let (g n ) be a bounded sequence in Y) and weakly compact 

in L 1 (}0, T[xB R xY R ), for allR> 0. Let (G n ) be a bounded sequence in L 1 (]0, T[xB R x 
Y R ) f or all R> 0. We assume that 

d t g n + —-Vx9n = G n in V(]0, +oo[xR 3 x Y). 
m 

Then, for any function \l/ G L°°(Y 2 ), with compact support, the sequence 



g n (t,x,y)^(y,y*)dy 

Y J n gN 

is strongly compact in L 1 (]0,T[xB R x Y R ), for all R > 0. 
This result can be improved: 

Corollary 3.10 With the assumptions of theorem \3.9\ we also have: 
for all R > and for any function * G L°°(]0, T[xB R x yj), the sequence 

g n (t,x,y)^(t,x,y,y*)dy 
is strongly compact in L 1 (]0,T[xB R x Y R ). 

Proof : The case of separated variables is obvious. Then, we proceed by a density 
argument as in |3j. 

□ 

Corollary 3.11 With the assumptions of theorem \ 3.9\ we also have: for all R > 
and for any sequence (^ n ) bounded in L°°(]0, T[xB R x Y R ) which converges a.e to 
G L°°(]0,T[x B R xY r ), the sequence 

g n (t,x,y)^ n (t,x,y)dy 

Y J n &\ 

is strongly compact in L 1 (]0, T[xB R ). 
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Proof : Let e > 0. The sequence (g n ) being weakly compact in L 1 (]0,T[x5^ x Yr), 
there exists S > such that 

WE G £(]0,T[x£fl x Y R ), \E\ < 5, sup / / / \g n \dtdxdy < e. 

n J J J e 

Then, by Egoroff theorem, there exists E G B(]0, T[xB R x Yjj) such that \E \ < S and 
\I/ n converge uniformly to \I/ on £7i := (]0, T[xBr x Yr) \ _E . Whence 



/ g n ^ n dy- / fi^A/ 
'y r Jy r 




< I I I \g n \\^n-^\dydxdt 

L 1 Q0,T[X-Bi{) ^0 Jb r Jy r 

< 2C e + sup - ^| /// \g n \dydxdt 

E-i J J J 



2C e+ oil) . 



We infer 



g n ^ n dy- I g n ^dy 

Yr JYr 



0. 



Li(]0,T[xB fl ) 



The sequence (^j g n ^/dy^j being compact in L 1 (]0, T[x5 fi ) in view of corollary 13.101 
the results follows. 

□ 

Now, we are able to establish the strong compactness of the sequence of f n ^-averages, 
and also the (Lf n ) one. 

Lemma 3.12 For all R > 0, and for all function ^ G L°°{Y) with compact support, 
f n (t,x,y)*(y)dy^ J f(t,x,y)*(y)dy m L x (\U,T[xBr). 

Proof : Since it is not clear that (Q~ (f n , f n )) is bounded in -^QO, T[xB R x Yr), we 
can not directly apply theorem 13.91 to the sequence (/„). For v > 0, we consider the 

sequence g v n :— — log(l + z//„) and we set 



qu , = Qt(fnJn) _ Q-UnJn) + <ff(/n) _ Qf (fn) 
1 + uf n l + vfn 1 + vfn 1 + vfn ' 

By the assumptions on (/„), we have 

d t g v n + ^-.V x g v n = G v n in T>'(]0, +oo[xIR 3 x Y). 



(3.50) 
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Since < g v n < f n , the weak compactness of (/„) established in the lemma [33] im- 
plies that (g^) is also weakly compact. Similarly, the sequence (g%) is bounded in 
L^jO^xM 3 x Y). Then, by corollaries ESI O and lemma EH (G v n ) is bounded in 
L 1 (]0, T[x5/j x Yr). Therefore, theorem 13.91 applies to (g£) for all v > 0. In particular, 
for all function \1/ G L°°{Y) with compact support and for all v > 0, the sequence 



g^(t, x,y)^(y)dy J is compact in L 1 (]0, T[xi?^), thus, by the uniqueness of weak 

Y Jn 

limit, 

J^g»{t,x,y)y{y)dy^ J g»{t,x,y)y{y)dy in L\]0, T[xB R ), (3.51) 

where g v is the weak limit of (g£) (up to an extraction). 
The result follows because 



sup sup / / \g v n - f n \dydx — -> 0, (3.52) 

n te[0,T] J JRSxY 

which implies the strong compactness in L 1 (]0, T[xi? fi ) of the sequence 

Uf n {t,x,y)^{y)dy^j . 

To show f)3.52p . we can use the inequality 

VM > 0, < u - i log(l + vu) = ^-ul {u < M } + ul {u>M} . (3.53) 
Then we obtain, for all n and for all t G [0,T], 
// \9n-fn\dydx < ^-f[ fndydx+ll f n i {fn>M }dydx 

We conclude by letting v — > 0, and M — > +oo. 

□ 

Proposition 3.13 We set p n (t,x) := J f n (t,x,y)dy and p(t,x) := J f(t,x,y)dy. 
Then, up to a subsequence, we have, for all R > 0. 

Pn — ► P in I 1 (]0,T[xB fl ) and a.e. (3.54) 

Proof : We have 

Pn = Pn+°n, wh ere pn ■= f n (t,x,y)dy. 

JYu 
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By the preceding lemma, p n 
M > 0, and 



M 



f(t,x,y)dy in L x (]0, T[xB R ) for all 



f n (t,x,y)dy < 



Y—Yil, 



< 



M 



Cte 



Y M 



y-y m 



Y 



\y\fn(t,x,y)dy 



2m 



m+ - he f n (t,x,y)dy, 



whence cr* f — > in L 1 (]0, T[x5 R ), uniformly in n. 



M—>+oo 



□ 



Lemma 3.14 For all R > 0, we have, up to a subsequence, 

Lf n — > Lf in LWO^IxBrxYr) and a.e. (3.55) 

Proof : Applying the corollary 13.101 with *$f(y,y*) = A(y,y*)i^ yE Y R }^{ y *eY R *}, we infer 
that the sequence (^j g^(t,x,y*)A(y,y*)dy*^j is compact in i^QO, T[xBr x Y r ) 
for all R* > 0. Using (13.52f) again, we obtain, for all R* > 0, the compactness of 
{^j fn{t,x,y*)A(y,y*)dy* S j in ^QO, T[xB R x Y R ). We conclude similarly as for 
the proof of lemma [3761 establishing 

lim sup / Af*dy*- / Af*dy* 

R*-++oc n J Yr ^ J y 

and identifying the weak limits. 



L l Q0,T[xB R xY R ) 



□ 



3.5 Regularity in time of the limit / 

In this subsection, we show the continuity in time of the limit /, which gives a sense 
to the Cauchy data /(0) = f°. 

Proposition 3.15 In fact, we have f G C([0, T],L 1 (R 3 x Y)). 

Proof : We use the integr al for mulation. Each g v n is a distributional solution of the 
renormalized equation, by (13.501) . so a mild solution. Therefore we have, for a.e (x, y) G 
R 3 x Y, and for all t, t + h G [0, T], 

/t+h 
G v n (a,x,y)da, 
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thus 

\K\t + h)-g»J(t)\\ L1{BRxYR) < !! f + \G»(*,x,y)\da. 

J JB R xY R Jt 

Moreover, by the subsections 13.21 and 13.31 the sequence {G u n ) n is weakly compact in 
L x {\Q,T[xB R x Y R ), thus for all t e [0,T], 

lim sup \\g$(t + h)- 9 u J(t)\\LHB R xY R ) = 0. 

Therefore the sequence (g$) is equicontinuous in C([0, T], ^(Br x Yr)). By the com- 
pactness of [0, T], this sequence is in fact uniformly equicontinuous, thus 

lim sup sup \\g v n \t + h) - ^(t)|Ui(s HX r«) = 0. 

ft->0 n ie [o,T] 

Then, (13321 and the estimate (I3T45D yield 

lim sup sup \\fl(t + h)-fl{t)\\ LHM .s xY) = 0. 

n te[0,T] 

Ascoli theorem entails that the sequence (/|) is compact in C([0, T], L 1 (IR 3 x Y)). 
The uniqueness of the limit in V'QO, +oo[xM 3 x Y) yields f G C([0, T], L\R 3 x Y)), 
and so / G C([0, T], L X (R 3 x 7)) by change of variables. 

□ 



3.6 Passing to the limit in a new integral equation 

Even if the renormalization provides weak compactness, a new problem appears: we 
will not be able to pass to the weak limit in (ECFR), because of the non-linearity of 
the factor f n / (1 + /„). That's why we need another formulation to our problem, which 
avoids the renormalization. But, remember that the term Q~(f, f) can not be defined 
as a distribution, so w e will use an integral equation which doesn't involve this term. 
We proceed as in [10] . 

We denote by T the linear transport operator 

T = d t + ^V x . 
m 

Let T -1 be the resolvant of transport operator, defined by: for g(t,x,y), we set u = 
T~ x g if u\ t=0 = and Tu = g. So, we have 

T~ l g(t, x, m,p, e) := / g(s,x — (t — s)p/m,m 1 p 1 e)ds. 
Jo 

T _1 satisfies the following properties: 
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(i) Forall R > 0, T^ 1 (L\]0, T[xB R x Y R )) C C([0, T], L X (B R x Yr)) continuously 
and weakly continuously. 

(ii) T~ l is nonnegative (Vg > 0,T _1 g > 0). 

Forall F e C([0,T], L l (B R x Yr)) such that TF > 0, we set 

Tp 1 = e~ F T~\ F . 

This operator is well defined from Z^QO, T[xB R x Y R ) to C([0, T], L^Sr x Yr)) and 
has the same continuity properties as T~ . 

Moreover, if (F n ) is a bounded sequence in C([0, T], L 1 (B R x Yr)) such that TF n > 0, if 
F n (t, x, y) — - > F(t, x, y) for all t and a.e (x, ?/), and if g n — ^ (7 weakly in -^ 1 (]0, T[x5 B x 
Yr), then 

Vte[0,T], Tp^g n (t) — Tp 1 g(t) weakly in L 1 ^ x Yr). 

The operator T^ 1 allows us to build an new formulation of our problem, which is better 
because it only involves Qt{f n ,fn), Q~f{fn), Qc(fn)- 

Lemma 3.16 / G C([0, T], /^(IR 3 x Y)) is a mild solution of (ECF) with initial data 
f(Q) = f° if and only if 

f = e- F f°(x - tp/m, y) + T/(Q c + (/, /)) + T F \Q + f (f)) - T/(Q 6 7(/)), (3.56) 
where F := T'^Lf). 

Proof : The result is deduced from the following fact: if / is a distributional solution 
of (ECF), then 

T(e F f) = TF e F f + e F Tf = e F (fLf + Q+(/, /) - Q~(f, f) + Q+(f) - Qj(/)) 

= e F (Qt(f,f) + Q J (f)-Q](f))- 

□ 

Now, we can finish the proof of theorem I3.2L 

End of the proof of theorem 13.21 We will pass to the weak limit in the following equa- 
tion, satisfied by each (/ n ): 

fn = e~ F "f°(x - tp/m,y) +T Fn \Qt(f n J n )) + T Fn \Q+(f n )) - T Fn \Q;(f n )), (3.57) 
where F n := T~ l (Lf n ) 

Notice that in view of (13.550 and the continuity properties of T~~ , the sequence (F n ) is 
bounded in C([0, T], L 1 (-B_r x Y r )), and F n (t,x,y) — > F(t,x,y) for all t and a.e(x, y). 
Thus, we can pass to the weak limit in the terms Tp n (Q~^ (f n )) and Tp n (Q~ (f n )) thanks 
to the lemma [3781 The term e~ Fn f®(x — tp/m,y) can be treated with the continuity 
in t = (for the L -norm) of each and /, established in the previous section. 
Eventually, the last term Tp (Q+ (f n , f n )) also pass to the weak limit thanks to the 
following lemma and proposition, which use the a.e convergence of the y-averages 
obtained in the previous subsection. 
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Lemma 3.17 For all R > and for all function ip G L°°(]0, T[xB R x Y r ), we have, 
up to a subsequence, 

Sy Qtifn, fn){t, x, y)(f{t, x, y)dy n J Y Q+(/, f)(t, x, y)<p(t, x, y)dy 
l + p n (t,x) l + p(t,x) 

in L 1 (]0,T[x.B fl ) and a.e. 
Proof : We have 

f Y Qt(fn,fn)<pdy 1 f x H f n (t,x,y)A(y,y*)<p(t,x,y + y*)dy \ ^ 



1 + Pn 2 J Y " ' Vy 1 + P„(t, S) 

Now, we apply the corollary 13.101 with ty(t,x,y,y*) = A(y,y*)ip(t,x,y + y*) 
(notice that \I> G £°°(]0, T[xB R x Y 2 2 K ) thanks to (Oil ). 

Therefore, the sequence g'^{t,x,y)A(y,y*)ip{t,x,y')dy S j is compact in 
L 1 (]0,T[x.Bfl x Yr), and we have 

x, x, y')dy / ^(t, x, y) y*)^(t, x, y')*/ (3-58) 

y </y 

in L 1 (]0 ) T[xS i i x Y R ), for all i/ > 0. 

Using (13.521) again, we obtain, up to a subsequence, that 

f n {t,x,y)A(y,y*)ip(t,x,y')dy / f(t,x,y)A(y,y*)<p(t,x,y')dy (3.59) 
y Jy 

in L\}0,T[xB R x Yr) and a.e. 

Up to another extraction, we infer, by (13.591) and (13. 54ft . 

fn(t,x,y)A(y,y*)ip(t,x,y')dy „ f f(t, x, y)A(y, y*)ip{t, x, y')dy 



Y 



l + p n (t,x) J Y l + p(t,x) 

a.e in (t,x,y) e]0,T[xB R x Y R . 

fn(t, x, y)A(y, y*)<p(t, x, y')dy 



Applying the corollary 13.111 with ty n (t,x,y*) :- , 

Jy 1 • pJ'-x) 

(which satisfies the required assumptions because ip is compact supported and A is 
locally bounded), we obtain the compactness of the sequence 



/ g v n {t, x, y*)V n (t, x, y*)dy* j in Z^QO, T[xB R ), and so, by f l332l . we deduce that 

/ f n {t, x, y*)^ n {t, x, y*)dy* J is compact. 
Jy J neN 

Finally, we conclude that 

J f n (t, x, y*)* n (t, x, y*)dy* ^ J y /(*, x, y*)*(t, x, y*)dy* in L\}0, T[xB r ). 
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□ 

Proposition 3.18 Up to a subsequence, we have, for all R > 0, 

QUfnJn)^Qt(fJ) weakly in L\]0, T[xB R x Y R ). 
Proof : We know by the lemma [3T4l that there exists Q(t, x, y) such that for all R > 0, 

Qi(fn,fn)^Q weakly in L^T^BrxYr). 

By (13.541) and a standard integration argument (we can refer to [17] for a proof), it 
leads to 

Qt } fnJn) ^-g- weakly in L\]0, T[xB R x Y R ). 

1 + Pn 1 + P 

Moreover, the previous lemma shows that 

Qtifnifn) Q~cifif) ii • rl/ln tT., d w N 

— — ^ — — weakly in L [\0,T[xB R x Y R ). 

1 + Pn 1 + P 

We conclude identifying weak limits. 

□ 

We have shown that / is a mild solution of (ECF). Since Qj(fn,fn), Qf(fn) and 
Qj(fn) converge weakly to Qt(f, f), Q~f (f) and Q](f) respectively, these three terms 
lie in L loc , and a fortiori, 

Qt(fJ) Q + fU) QJ(f) r r i 



i + / ' i + /' i + / 



eL; oc (]0,+oo[xR d xF) 



The term — | — '— is automatically in L; oc because Lf G L (]0, T[x.Br x Yr) for all 

i2 > o and 9lillll < L /. 

1 + / " 

Thus, / is indeed a renormalized solution of (ECF). 

□ 
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